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Abstract 



In this paper, we consider a moving rigid solid immersed in a potential fluid. The fluid-solid 
system fills the whole two dimensional space and the fluid is assumed to be at rest at infinity. 
Our aim is to study the inverse problem, initially introduced in [ '], that consists in recovering 
the position and the velocity of the solid assuming that the potential function is known at a 
given time. 

We show that this problem is in general ill-posed by providing counterexamples for which 
the same potential corresponds to different positions and velocities of a same solid. However, 
it is also possible to find solids having a specific shape, like ellipses for instance, for which the 
problem of detection admits a unique solution. 

Using complex analysis, we prove that the well-posedness of the inverse problem is equivalent 

to the solvability of an infinite set of nonlinear equations. This result allows us to show that 

. . when the solid enjoys some symmetry properties, it can be partially detected. Besides, for any 

solid, the velocity can always be recovered when both the potential function and the position 

k> are supposed to be known. 

%_j Finally, we prove that by performing continuous measurements of the fluid potential over a 

j^ time interval, we can always track the position of the solid. 
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1 Introduction 

1.1 History 

Sonars are the most common device used to spot immersed bodies, like submarines or banks of 
fish. These systems use acoustic waves: active sonars emit acoustic waves (making themselves 
detectable), while passive sonars only listen (and therefore are only able to detect targets that are 
noisy enough). To overcome these limitations, it would be interesting to design systems imitating 
the lateral line systems of fish, a sense organ they use to detect movement and vibration in the 
surrounding water. 

Most of the published results on inverse problems in Fluid Mechanics concern the detection of 
fixed immersed obstacles. For example in [ ] the authors prove that a fixed smooth convex obstacle 
surrounded by a fluid governed by the Navier-Stokes equations can be identified via a localized 
boundary measurement of the velocity of the fluid and the Cauchy forces. In [ ], the authors 
identify a single rigid obstacle immersed in a Navier-Stokes fluid by measuring both the gradient of 
the pressure and the velocity of the fluid on one part of the boundary. The distance from a given 
point to an obstacle is estimated in [ ] from boundary measurements for a fluid governed by the 
stationary Stokes equations. 

To our knowledge, the only work addressing the detection of moving bodies is [3]. In this 
paper, the authors consider a single moving disk in an ideal fluid and prove that the position and 
velocity of the body can be deduced from one single measurement of the potential along some part 
of the exterior boundary of the fluid. They obtain linear stability results as well, by using shape 
differentiation techniques. 

1.2 Problem settings 

Domains, frames, coordinates 

At a given time t, we assume that a rigid solid occupies the domain S C R 2 , while the domain 
J- := R 2 \ S is filled by a perfect fluid. Let us assume that S is a simply connected compact set. 
The unitary normal to dT directed towards the exterior of T is denoted by n. As being a rigid 
solid, S is the image by a rotation and a translation of a given reference domain 5 which will 
merely be called in the sequel the shape of the solid. Therefore, at any time, there exist an angle 
8 e R/27T, a rotation matrix R(9) £ SO(2) of angle 8, a point s := (si,S2) T € R 2 (the center 
of the rotation) and a vector r := (ri,r2) T of R 2 such that S = R(9)(Sq — s) + r. We will be 
concerned with recovering the position of the solid, so it is worth remarking that the triplet (8, s, r) 
is not unique. Two triplets (8j,Sj,rj) E R/2ir x R 2 x R 2 (j — 1,2) give the same position for 
any iSo if and only if R(8i) = i?(#2) = R and R{s\ — 82) = ri — r2. These equalities define an 
equivalence relation in R/27T x R 2 x R 2 . However, we want also to take into account the possible 
symmetries of the solid. So, given So, we say that two triplets (8j,Sj,Tj) are equivalent when 
R($i)(So — si) + n = R(02){So — s 2) + i"2- We denote by V the set of all of the equivalence class 
p. We will make no difference in the notation between p and any element (8, s, r) belonging to this 
class. In particular, we will write in short that for any x £ R 2 , px = R(8)(x — s) + r. In the sequel, 
p will be merely referred to as position of the solid. 

Later on, we will use tools of complex analysis, so rather than R 2 , we will sometimes identify 
the plane with the complex field C. For any complex number z := z\ + 1Z2 (i 2 — — 1, Zi,z% € R), 
we will denote z :— z± — iz2 the conjugate of z and D will stand for the unitary disk of C. 



Sequences of complex numbers 

For any sequence of complex numbers c := (ck)kez, we can define c := (ck)kez and c := (c-k)kez- 
For any two sequences a := (a^kez an d b '■= (&fc)fceZj we recall the definition of the convolution 
product: a*b := (X^ez a k-jbj)kez- The convolution product can be iterated n times (n an integer) 
to obtain a n := a * a * ... * a. 

Rigid velocity 

The solid is moving. We denote by v(x) := (vi(x), V2(x)) T G R 2 the rigid Eulerian velocity field 
defined for all x G S. This notation turns out to be v(z) :~ «i(z) + iv^z) in complex notation. 
It is well known in Solid Mechanics that v can be decomposed into the sum of an instantaneous 
rotational velocity field and a translational velocity field. For any x :— (x\, X2) T G R 2 , we introduce 
the notation x := (— X2, xi) T and we have v(x) = lj(x — s) 1 - + w, where s <G R 2 is the center of the 
instantaneous rotation, u> £ R is the angular velocity and w := (wi,W2) T G R 2 the translational 
velocity. Since we wish to recover these data, it is worth observing that the triplet (w, s,w) G 
R x R 2 x R 2 is not unique. Both triplets (ujj, Sj, Wj) (j = 1, 2) give the same rigid velocity field v 
if and only if lu\ = u>2 = to and w(s^ — s^) = wi — W2 (in particular, we can always choose for s 
any point of R 2 ). This is an equivalence relation and the velocity v can be seen as an equivalence 
class. We denote V as the set of all of the equivalence class and we will not differentiate, in what 
follows, between the vector field, the class of equivalence and any element of this class. All of them 
will be denoted by v. 

Definition 1.1 (Configurations). For any given shape Sq, we define a configuration as any position- 
velocity pair (p,v) ePxV. 

Fluid dynamics 

The dynamics of the fluid is described by means of its Eulerian velocity field u(x) :— (ui(x),U2(x)) T 
defined for all ieJ. Since the fluid is assumed to be perfect (i.e. incompressible and inviscid) and 
the flow irrotational, there exists a potential function ip, harmonic in J 7 , such that u(x) = V<^(a;) 
(x G J-). The fluid is assumed to be at rest at infinity so we impose the asymptotic behavior 
|V(/3(x)| — S- as \x\ — *■ +00. The classical slip boundary condition for inviscid fluid reads as 
u • n = v • n on dS and yields a Neumann boundary condition for ip, namely d n ip = v • n on dS. 
Although the domain T is not simply connected, we can still consider tp, the harmonic conjugate 
function to tp, because J ds d n pda = 0. The functions ip and ip satisfy the relation Vi[> = (Vip) 1 - in 
J- . In Fluid Mechanics, -0 is called the stream function and the complex function £ = ip + iip is the 
holomorphic complex potential. As usual, we define u := u\ +iu2 = £' as the complex fluid velocity. 
Observe that the complex potential, as being solution of a boundary value problem, depends on the 
domain T and the velocity v only. With the notation introduced earlier, we deduce that £ depends 
only on the shape S and the configuration (p,v) ePxV. 

The complex potential is defined up to an additive constant which can be chosen such that 
|£(z)| —¥ as \z\ — >• +00. For any i/£C, the complex potential can be expanded in the form of a 
Laurent series: 



^ (z - v)i ' 



£(*)== 2.7^7' \z-v\>R{v), (1.1) 



where Xj(f) (j > 1) are complex numbers and R{v) '■= limsup ? _ i . +00 | Xj (v) | 1 ' i '. The series is 
uniformly convergent on the set {z G C : \z — v\ > R(is)}. 

Measurements 

We measure the complex velocity u of the fluid in some open subset of J- ' . The Analytic Continuation 
theorem tells us that we can deduce the value of £' everywhere in the connected open set T and 
then also the value of £, up to an additive constant. In particular, we will assume that for all v G C, 
we can always evaluate all of the terms of the complex sequence (A(f))j>i arising in the expression 
(1-1)- 

1.3 Main results 

Definition 1.2 (Dctcctability). A solid of shape Sq is said to be detectable if, for any configuration 
(Pi v ) £ V x V, the knowledge of the potential holomorphic function £ suffices for recovering the 
pair (p,v). 

Observe that this definition makes the property of being detectable independent of the configu- 
ration: detectability is a purely geometric property of the solid. Our first result is that not all the 
solids are detectable: 

Theorem 1.3. For any integer n > 2, there exists a holomorphic function £, a shape Sq, and 
n configurations (pj,Vj) G V X V, j = 1, . . . ,n satisfying pj 7^ pt if j 7^ k such that £ is the 
potential of the fluid corresponding to the solid of shape Sq with any of the configurations (p,-, Vj), 
j = l,...,n. 

In other words, for any integer n, there exists at least one solid that can occupy n different 
positions with n different velocities and for which the fluid potential is the same. This theorem 
shows that the result obtained in [3] for a disk can not be generalized to any solid. However, not 
only the disk is a detectable body: 

Proposition 1.1. Any ellipse is a detectable solid. 

Going back to the general case, it is easy to see that the holomorphic potential never admits 
an analytic continuation over the whole complex plane. Furthermore, for any analytic continuation 
of the potential inside the solid, we will prove that the location of the singularities provides clues 
allowing one in many cases to determine the position of the solid. This discussion is carried out in 
Subection 6.1. 

According to Theorem 1.3, the problem of detection is ill-posed in the general case. However, we 
claim that when the solid enjoys some symmetry properties, it can be partially detected (i.e. some 
but not all of the parameters among r, a, w,cj can be deduced from the potential). The following 
proposition illustrates this idea: 

Proposition 1.2. If the shape of the solid is invariant under a rotation of angle tt/2 then r, w 
and \u>\ can be deduced from the potential function. 

We refer to Propositions 6.4 and 6.5 for a more precise statement of this result. 
In the general case, we can also try to determine less parameters with more information. For 
instance, we can prove: 



Proposition 1.3. For any solid with configuration (p, v) G V x V, the knowledge of both the 
potential function and the position p suffices for recovering v. 

Finally, we can also measure the potential function, not only at a given instant, but over a time 
interval. In this case, we obtain: 

Theorem 1.4 (Tracking). For any solid So, if we know its position at the time t — and we 
perform continuous measurements of the complex potential over the time interval [0, T] for some 
T > then we can deduce the configuration of the solid at any time t £ [0, T]. 

1.4 Outline of the paper 

In the next section, we provide examples of non-detectable solids and prove Theorem 1.3. In 
Section 3, we derive the expression of the complex potential. In Section 4, we determine all the 
stealth solids, i.e. all the solids that can move in the fluid without disturbing it. The detection of a 
moving ellipse is discussed in Section 5. Section 6 is split into three parts: the first one is dedicated 
to the study of the singularities of the potential function and the second one to its asymptotic 
expansion and how these results can be used for the detection problem we are dealing with. The 
third subsection deals with an example of detection. In Section 7 we give the proof of Theorem 1.4 
and at last in Section 8, we indicate some remaining open problems. 

2 Examples of Non-detectable Solids 

This Section is mostly devoted to the proof of Theorem 1.3. 

Expression of the stream function 

Let a shape Sq and a configuration (p, v) be given with v = tu(x — s) + w (for some real number 
uj and some vector s) and remember that S — p(<So). Then, let us introduce 7 : [0,£[n- -f(s) = 
(7i(s),7 2 (s)) T € R 2 a parameterization of OS satisfying |7'(s)| = 1 for all s € [0,£[ (£ > 0). We 
assume that dS is described positively (counterclockwise parameterization), we denote t = 7' (the 
unitary tangent vector to dS) and we get n = t- 1 . We deduce that d n tp = —d T i[> and hence that 
9 T 1/7(7) = — wi72+W27i+w7'-(7 — s). We can integrate along dS to obtain ^(7) = — W172+W271 + 
(w/2)|7 — s| 2 + C on dS, where C is real constant. This Dirichlet boundary condition for the stream 
function reads also: tp(x) = —w\X2 + W2^i + ( w /2)|x — s| 2 + C on dS. In this form, the boundary 
of the solid turns out to be a level set of the function g(x) :— (u/2)\x — s| 2 — W1X2 + W2X1 — ip(x), 
an observation we will now take advantage of. 

Proof of Theorem 1.3 

Pick some integer n > 2 and consider the harmonic function whose expression in polar coordinates 
is ip(r, 9) := cos(n9)r~ n . Since, in Cartesian coordinates, \(difj/dxj)(x)\ < \Vifj(x)\ = n|x| - " -1 (j — 
1, 2), we deduce that for any u > and s £ R 2 , there exists 6 > such that (dip/dxi)(x)— U)(xi — Si) 
and (dip / dx2){x) — lo(x2 — 82) can not be simultaneously null providing \x — s| > S. Applying the 
Local Inversion Theorem, we deduce that for any A € R, the solutions of 

~\x-s\ 2 -i>(x)-\ = 0, (2.1) 



satisfying \x — s\ > 6 (if any) are locally smooth curves. 

From the estimate 1^(^)1 < M™ ( x <= R- 2 ): we deduce that for all s := (si,S2) T and all 
u) > e > 0, there exists 5' > such that: 



A < -\x- 

- 2 1 



•V(«) 



\ ^ w + e i 
A < — - — kc 



for all A € R providing |x — s| > 5'. If we choose for instance A > max(<5', S) 2 (ui + e) 
level set of the function g(x) := u>\x — s| 2 /2 — ip(x) — A between the circles 



there is a zero 

= Vzx/V^^ 

and \x — s| = ^/2\/y/u) + e (because s/2\/y/uT+s > \/2S > S). It remains to choose s properly, 
in order to take advantage of the symmetry of the function ip. Let p be any positive number and 
denote S\ the zero level set of g obtained as described above by specifying si := (p, 0). This level 
set defines the smooth boundary of a solid for which ip is the stream function (and £(z) := i/z n 
the holomorphic potential) associated with the velocity vi := uj(x — Si) . By choosing next 
s fc = (pcos(2(k — l)7r/n), psin(2(fc — l)7r/n)) T for k = 2, . . . ,n, we obtained n — 1 copies of Si at 
n — 1 different positions with respective velocities v^ := oj{x — s&) . Some examples of such solids 
with the associated rigid velocity fields are displayed in Figures 1, 2 and 3. 





(a) w = 0.7, si = (0.7,0), A = 1.5 



(b) w = 0.7, s 2 = (-0.7,0), A = 1.5 



Figure 1: For both configurations, the stream function is the same. It reads ip(r,9) := cos(2#)/r 2 
in polar coordinates. The holomorphic potential is £(z) := i/z 2 . 



3 The Complex Potential 

Before going further, we need to describe the shape Sq. Actually, for convenience, rather than Sq 
we shall describe Tq := C\S~q. Thus, assume that .Fo is the image by a conformal mapping / of 
Q := C\D, the exterior of the unitary disk. For any simply connected shape Sq and corresponding 
domain J-q, the Riemann Mapping Theorem tells us that / can be written in the form: 



E 
fc<-i 



f(z) = ciz + cq + V c k z k , (ze£l), 



(3.1) 



where c k & C for k = 1 and all k < —1 and C\ ^ 0. We can assume, without loss of generality, that 
Co = 0. To simplify forthcoming computations, we will also assume that c k is actually defined for 
all k € Z and that c k = for k = and k > 2. We denote c := (ck)k^z the complex sequence of 
elements c k and the Area Theorem (see [ , Theorem 14.13]) tells us that the area of Sq is equal to 
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-3-2-10123 
(a) u) = 0.6, p = 0.5, A = 1.5 and si = 
(pcos(7r/4), psin(7r/4)) 




-2-1 1 2 

(b) w = 0.6, p = 0.5, A = 1.5 and s 2 

(pcos(— 7r/4),psin(— 7r/4)) 



Figure 2: For both configurations, the stream function and the holomorphic potential are the same 
as in Figure 1. 



7r Sfc<i ^l c fe| 2 - Since Sq is of finite extent, it means that this sum has to be finite. Actually, we will 
assume also that c € £ (C), which entails in particular that / is continuous in the closed set Cl. 

Such a description allows us to consider a broad set of solids. In particular, the boundary of 
the solid can be very rough. Degenerate cases can be considered as well (for instance Sq can be a 
segment modeling a one dimensional beam). 

For any position p := (R(a), 0, r), we recall that S :— R(a)So + r is the actual domain occupied 
by the solid. Let us introduce then the functions ipo(x) := tp(R(a)x + r) and ipo{x) := ip(R(a)x + r) 
which are harmonic (and defined) over the fixed domain J-q, For any velocity v := (w,r, w) (we 
choose here s = r), the Dirichlet boundary condition for -0 turns out to be, in complex notation 
2i?/> : = Woz — wqz + iuj\z\ 2 where Wq(z) := w(R(a)z + r). We introduce £ the holomorphic complex 
potential of the fluid defined for any z £ il by ((z) — ipo(f(z)) + ii/jo(f(z)). Since z — \j z on dil, 
we get the identity 2iipo(f(z)) = —wf(z) + wf(l/z) + icuf(z)f(l/z). For any z £ dfl, we have also 
/(!/*) = E/cgz c-kz k = Efcez CkZ k and f(z)f(l/z) = Efeez( a * c) k z k . So we get: 



2iip (f(z)) = ^[-w Cfe + w a c k +iuj(c*c) k ]z k , (z G dtt). 



(3.2) 



fcez 



According to [ 6, Chap. IX, §9.63], we keep only the negative powers in (3.2) to get the expression 
of C- Defining the coefficients Ck{wo, w) := [— WqCu + w$c k + iu>(c * c)k] for all fc < — 1, we obtain: 



CC0= E CfcK,^)z fe , (zefl). 



fe<-i 

Eventually, the expression of the measured complex potential, defined in J-, is: 

az) = ar 1 ((z-r)e-n), (z&n 

According to our rule of notation, we introduce as well 

tow = <po(z) + Mz) = a.r\z)), ( z e j- ). 



(3.3) 

(3.4) 
(3.5) 
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Figure 3: The stream function is ijj(r,6) = cos(6#)/r 6 , the holomorphic potential is £ = i/z 6 and 
cj = 0.7, p = 0.9, A = —2.5 and Si = (pcos(fc7r/6),psin(2fc7r/6)) for k = 1, . . . ,6. 



4 Stealth Rigid Solids 

In this Section we wish to determine all the possible shapes and configurations of solids for which 
the complex potential £ is identically null. Such a displacement will be termed stealth. 

Theorem 4.1. The only solids S that can undergo stealth motions in a fluid are: 

• Disks rotating about their centers; 

• Arc of circles and segments with velocity field everywhere tangent to S. 

The arc of circles and segments are one dimensional solids and can be considered as degenerated 
cases. 

Proof. Let us assume that £ = 0. Then we have also £ = which means that —2^s(uf(z)) + 
uj\f(z)\ 2 = for all z € dft. If u ^ 0, some easy computations tell us that for all z € dCl, f(z) 
belongs to the circle of center iwo/2uj and radius |iuo|/(2|u;|). Since / is an homeomorphism from 
dil onto f(dfl), f(dil) is a connected compact subset of this circle. 



• If f(dfl) is the complete circle, it means that c\ = 1 and Ck = for all k ^ 1. In this case, 
since £i = and (c* c)i =0, we deduce that wo — and hence that the circle is just rotating 
about its center. 

• Up to a translation and a rotation, all the conformal mappings that map the circle onto an 
arc of circle have the form f(z) = z + (1 — h 2 )/(z + ih) where h is any real number such that 
< h < 1. We can put / into the general form (3.1) by setting: c^ — 1, c_i = 1 — h 2 and 
Ck = (1 — h 2 )(—ih)~ k ~ 1 for all A; < —2. Some simple computations lead to: (c*c)_i = —ihc-i 
and (c * c)fe = (ih^ 1 — ih)ck for all fc < —2. Substituting these expressions into (3.3) and 
writing that £fc(u>o,w) = for all k < — 1 we obtain the same equation for all k which yields 
the relation: wq — uj{h — 1/h). We can then easy prove that this motion corresponds to the 
case where the velocity field is tangent to the solid. 

Let us assume now that uj is zero (and wq ^ 0). In this case, we deduce with (3.3), that Ck = for 
all k < —2. For k — — 1, we get c_i = ciWq/wq. We set wq — Re 10 , C\ — Re 1 ^ and we rewrite / in 
the form: 



/(*) 



R e i0 z 



,i(-0+26) 



I* 



2Re 



it) 



AP-8). 



-i(ff-e) 



I* 



We seek the image of the unitary circle by /. We specify z = e lt with t G R/2tt and we get 
f(e lt ) — 2Re 10 cos(/3 — 9 + t). So the image of the unitary circle is the segment [-R, R] turned by 
an angle 9. The velocity wq is collinear to the segment. D 



5 Detection of a Moving Ellipse 

When So is an ellipse, the function / has the form f(z) = (a + b)z/2 + (a — b)/2z, where a, b € R+, 
a > b > 0. We can now give the proof of Proposition 1.1. 



Proof. First, we can explicitly compute the inverse function 



_ {a 2 - b 2 ) 



(z e T Q ). 



(5.1) 



In this expression, —y/a 2 — b 2 and V 'a 2 — b 2 are branch points and the function is holomorphic 
everywhere but on the segment [-ya 2 — b 2 , %/a 2 — b 2 ] which is a branch cut. Next, we get: 



C(z) 



a — b 



w - 



w - 



a + b 



and then: 



£(*) 



[-(a 2 - b 2 )w + (a + 6) 2 w ]e i 



(zeQ), 



i(a 2 -b 2 )(a + b) 2 e 2la u 



2{z- 




(a 2 - b 2 )e 2 



(z-r) 2 



4(z - r)< 




1 - 



{a 2 - b 2 )e 
{z-rf 



(z e F). 



Observe that, due to the symmetry of the ellipse, we can change a into a + ir and accordingly 
wo into —wq without changing the expression of £. The potential £ is holomorphic everywhere 
but on the branch cut [r — \fa 2 — b 2 e la , r + v 'a 2 — b 2 e ta ]. So if we know thoroughly £, we can 
determine the location of the branch points r — yfa 2 — b 2 e la and r + \J a 2 — b 2 e %a and hence also 



the position of the center r and the orientation a (up to 7r only). Wc compute next the limit 
/i := lini|£|_ s . +00 e~ la t;(z)z/(a + b) = [—(a — b)w + (a + b)wo]/i and we deduce the expression of 
wo, namely: wq — (fi + p.)/b + (/i — p,)/a. The only remaining unknown quantity u is next easily 
obtained, following the same idea. □ 




Figure 4: Level sets of the holomorphic potential £oj for a = 2, b = 1, wq = e I7r / 3 and uj = —2. The 
boundary of the ellipse (dashed line) can hardly be directly detected but the branch points — \/3 
and \/3 are clearly identifiable. 



6 Detection: General Case 

6.1 Singularities of the holomorphic potential 

In the preceding example, the branch points of the potential £ played a crucial role in determining 
the position of the solid in the fluid. Notice that the existence of these points did not depend on the 
configuration but only on the shape of the solid (they came from the definition (5.1) of the inverse 
function / _1 and were subsequently just translated and rotated according to the position). We 
shall prove that this result can be generalized to any solid: there is no singularity in the potential 
function, that does not come from the conformal mapping / _1 (but unfortunately the potential 
function may have less singular points than the function Z" 1 ). Let us make this statement precise: 

Definition 1 (Analytic continuation). An holomorphic function £ (respect. £oJ defined in a con- 
nected open set JF (respect. J-q) containing T (respect. J-q) is called an analytic continuation of £ 



(respect. S,q) when £, — S, in T (respect. £o = £o in Tq). 



There may exist several analytic continuations of £ that do not coincide everywhere. Assume that 
£i and £2 are two such functions defined respectively on T\ and T-2- So the Analytic Continuation 
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Theorem ensures only that £1 = £2 on the connected component of F\ n F 2 containing F . In 
Section 5 for instance, we can not choose where the branch points are, but there are many different 
possible choices for the branch cut, each one corresponding to a different analytic continuation of 

e 

Assume that for some potential function £, there exists an analytic continuation £ such that 
F = C. Since, by construction, £(z) tends to as \z\ goes to infinity, £ is a bounded entire 
function. According to Liouville's Theorem, this function is constant, equal to 0. This case was 
treated in Section 4 and is possible only for solids listed in Theorem 4.1. For all of the other solids 
and for any analytic continuation £, there exists at least one point, located inside the solid, which 
does not belong to F. This very simple observation allows one to locate the solid in a very first 
approximation. 

Let us now prove that the singularities of £ come from the singularities of / . 

Proposition 6.1. If there exists an analytic continuation g of f^ 1 defined on an open connected set 
Tq containing Fq, then for any configuration (p, v) 6 V x V, there exists an analytic continuation 
f oft defined on F := p(F \ 3 _1 ({0}))- 

In this proposition, the notation g _1 ({0}) stands for the preimage of {0} under g and does not 
mean that g is invertible. Since g is holomorphic, the set 3 _1 ({0}) consists only in isolated points 
and Fq \ g~ 1 {{0}) is still connected and still contains Fq. 

Proof. For all z € dD, we can rewrite (, in the form: 

C(z) = -w (f(z) - C\z) + wqCiz' 1 + ioj c 1 z~ 1 (f(z) - c\z)+ 

J2c- k z k (f(z)- J2 <***)] ■ (6- 1 ) 

fc>l -k<j<l 

Expanding the right hand side and recombining terms, we get the identity: 

C(z) = -Wof{z) + ciWqz + uciz^ 1 + iu f(z)f(z" 1 ) - H\(z) , (z E dD), 

where Hi(z) := Ylj>o(c * c )j^ an< ^ f( z ) := ^i z + Tlk<-iCkZ k ■ Classical results for the con- 
volution product ensure that this series is uniformly convergent for \z\ < 1 since ||c * c||£i(c) < 

l|c||£i(c)||c||£i(c)- We next obtain that: 

C(/ -1 (*)) = ~w z + ^wof-^z) + uc 1 /f- 1 (z) + iu[zf0./f~ 1 (z)) - H.if-^z))] , (z G 8S Q ). 

Let g by any analytic continuation of / _1 (not necessary invertible), defined in an open set Fq. It 
can be split into three parts: Fq :— {z e Fq : \g(z)\ > 1}, Fq := {z e Fq : < \g(z)\ < 1} and 
Fq := {z e C : g(z) = 1}. We can next define: 

io(z) := - w z + ciw g(z) + w c\/g{z) + ioj zf{l/g{z))-Hx{g{z)) , zeFqUFq, 



Uz) 



w Q z + cxwogiz) + w ci/g(z) + iu H 2 {g(z)) 



z e F + , 



where H 2 {z) := Ylj<-i(c* c )o z ^ ^ s uniformly convergent for \z\ > 1. We deduce that the function 
£0 is holomorphic in Fq and in Fq and continuous in Fq \ <? _1 ({0}). Let z € Fq and denote 



11 



/ ^ 



m Z ) = o}- 



V" 




D{z ,R) 

Figure 5: The conformal mapping g around the point z for n = 3. 



z\ := g(zo). Since g is holomorphic at the point zq, there exists R > 0, n > 1 and a function h 
holomorphic in the disk D(0,R) such that h(0) ^ and 

g(z) =zi + (z- z ) n h(z - zo), z G D(0, i?). 

This identity allows us to describe the set T$ nearby the point zq. For instance, if n — 3, and for 
R small enough, we get something like in Figure 5 where z\ = g{z a ), U| =1 A+ = D(z ,R) fl Jq , 
U^ = i^47 = ^( z o^) H J^(7 and the curves radiating from z correspond to the set D(z ,R) flJj. 
Except maybe at the point z , the boundaries shared by the regions A^ and A~ are smooth. The 
function g maps A^~ onto C/ + (an open set located outside the unitary disk) and A~{ onto U~ (an 
open set located inside the unitary disk). Furthermore, the function ffl^+y^- : AjLiAj —t U~UU + 
is a conformal mapping. Consider next a conformal mapping cj) which maps a neighbourhood of 
z\ onto a neighbourhood of as in Figure 5 and such that the image of the unitary circle is the 
imaginary axis. We can then apply [ , Theorem 16.8]: the function £o ° 9 1 ° 4> _1 i s holomorphic 
on both sides of the imaginary axis and continuous across this boundary, so it is holomorphic over 
the whole domain. We deduce that £o is holomorphic across the boundary between AJ and A^ . 
We can repeat this process with the domains A^ U A^ , A\ U A^ and so on. Finally, we obtain that 
£o is holomorphic on the whole disk D(zq,R) except maybe at the point Zq. But once more, the 
continuity of the function £o does not allow this possibility. □ 



As already mentioned, the potential function can make some singularities of / to vanish. This 
is illustrated by the example of Figure 6. 

6.2 Asymptotic expansion of the holomorphic potential 

In this section, we shall compute the asymptotic expansion of £ in terms of the geometrical data of 
So and the configuration (p, v) G V x V. As explained in the preceding section, we know that if Sq 
is not one of the solid listed in Theorem 4.1, the potential function admits no analytic continuation 
on the whole complex plane. It allows one to deduce approximately where the solid is. For all 
v G C, we can next consider T, a contour large enough to encircle the solid and the point v. The 
contour T := f(e~ la (T — r)) encircles the unitary disk. According to the expression (1.1) of the 
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Figure 6: The solid consists in a disk of center C with a segment [A, B]. When this solid is moving to 
the left, parallel to the segment, one can easily check that the potential coincides with the potential 
of the disk. Although any analytic continuation of f^ 1 has singularities at A and B, the complex 
potential does not see these points. In this case, the singularities do not allow one to determine the 
orientation of the solid. 



potential as a Laurent series, we obtain that, for all n > 1: 

A» = ^- I £(*)(* - vT-Wz = ^1 C(z) (e ia f(z) + r - v)"' 1 f(z)dz 



1 1 

n 2i7r 



£ (f) e tka (r - v) n - k j>_ C(z)f(z) k dz. 



If we define the complex sequence d := (dk)k£Z by dk ■= {k + l)Cfc+i for all k < — 1 and dk = for 
k > 0, we obtain: 

a «w = -I E (l) etka ( r - »y i ~ k ( d * cfc )-i> (» e c )- 

We can rewrite the last term: 



n ' — ' \k 

k=l 



(d*c k u= ]T 



c ii ■ ■ ■ c i k di k+1 



ii + ... + i k + 1 =-l 

= -A k w + B k w Q + iuCk, 



where 



il+...+ifc + l=0 

C k := ^2 (*1 +*2)C-*iC« 



-Dfe- / ?lC_ i 1 Ci 2 . . . Ci k+1 , 



H+...+i k+1 =0 



2 ' ' • "-Jfc + 2 - 



il + ..+i fe+ 2=0 
i 1+ ; 2 <-l 



In the following, to simplify the notation, we consider the quantities Ak '■= —Ak/k, Bk '■— —Bk/k 
and Ck := —Cf./k. Indeed, we get, for all n > 1: 



fc=i 



n — 1 



a «m=e;:;^(-^ 



A»o + -Bfe^o + iwCfc 



(6.2) 
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The problem of detection can now be reformulated as a purely algebraic problem: the complex 
sequence (Aj(V))j>i being given for all v <G C, as well as the complex numbers Ak, Bk and Ck 
(k > 1), can we solve the infinite nonlinear system of equations (6.2) and find the values of r, a, Wq 
and w? According to the results of both Section 2 and Section 4, we already know that there exist 
cases (namely, coefficients (ck)k^z) for which the answer is negative. 

In order to rewrite this infinite set of equations in a convenient short form, we introduce some 
linear operators: let us denote by N any positive integer and define Gn ■ R 3 — ► C N by (GNU) k := 
i-A k + B k )Ux + i(A k + B k )U 2 + iC k U 3 for all U = (Ui,U 2 , U 2 ) T £ R 3 and all 1 < k < N. We 
define D N : C N ^ C N and S N : C N -> C N as well, by respectively (D N Z) k := kZ k (1 < k < N) 
and (S N Z) 1 := and (S N Z) k := Z k _ x (2 < k < N) for all Z := (Z u . . . ,Z N ) € C N . The first N 
equations (6.2) can now be rewritten as: 

A N (v) = e log(r-v)D N( ,S N D Ne - log(r-u)D Ne i a D N g NU ^ (g 3^ 

= @ N (r-u,a)g N U, (6.3b) 

where U := (K(w ),3(w ), w) T and A N (v) := (\i(v), . . . , \ N (is)) T . The operator e SjvDjv is lower 
triangular and the identity (e SNDN )k, n = (tZi) for all 1 < fc < n < N, not so obvious, can be 
found in [2]. Considering the expressions (6.3), it is worth noting that: 

• In (6.3), the coefficients Xj{v) in the asymptotic expansion of £ are obtained by applying to 
the vector U (the velocity) first the operator Gn encapsulating the information relating to the 
geometry of the solid and next the operator 0./v(r — i>, a) depending only on the position. 

• The linear operator Gn depends on the complex sequence c only, i.e. on the shape of the 
solid. Moreover, the complex quantities Ak — C-kc\ {k > 1), B k+ \ — c_ k cic\ (k > 2) and 
Cfc-i — c_feC__ic 1 "~ 1 (k > 1) do not depend on c_„ for all n> k. In other words, for all N > 1, 
Gn depends on c\, C_i, c_2, • • ■ , c_at_i only. We deduce: 

Proposition 6.2. Let Sq and Sq be two shapes described by means of the complex sequences 
(c\)k>\ and {c\)k>\, such that c\ = c\ for all 1 < k < N. Then, if both solids have the 
same configuration, their complex potentials will have the same asymptotic expansion up to 
the order N — 1 . 

• The solutions (r\,r2, a,^R.(wo), ( i(wo), ui) T of all of the equations (6.3) (for all N > 1), form 
a sub-analytic set of R 6 (because &n is analytic in r\ , r 2 and a) . This subanalytic set has 
dimension d with < d < 6. However, because the dependence in (#t(wo),^s(wo),uj) T is 
linear, if it had dimension d > 4, it would entail the existence of a position (ri,r 2 ,a) and 
a non-zero velocity U € R 3 such that 0Ar(r — v,oc)GnUq = for all N > 1. As already 
mentioned before, this case is only possible if the solid is in the list of Theorem 4.1. 

We do not know if there exist solids such that < d < 3. Observe that in Section 2, we have 
only given examples for which the solids can occupy a finite number of different positions, so 
d = in these cases. 

For all N > 1, we can invert the system (6.3) to obtain: 

GnU = e - iaDN e los{r - u)DN e- SNDN e- los{r -" )DN A N (v), (6.4) 

= G N (r-v,a)~ 1 A N (v), (6.5) 
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or equivalently, with the notation of equation (6.2), 



A n w + B n w + ilxlC-n 



for all n G N. In this form, we can easily prove: 



E 

,fe=i 



n- 1 
fc-1 



: X k (u) 



(6.6) 



Proposition 6.3. If the solid does not occur in Theorem 4-1 and its position is given, then we can 
deduce its velocity. 

Proof. Denote Qi (j — 1, 2,3) the complex sequences respectively defined by Q\ := (—A k +Bk)k>i> 
Q k := (i(Ak + £>fc))fc>i an d Q\ := (*Cfe)fc>i. If these sequences were not R-linearly independent in 
C N , it would exist (Uj_, E/ 2 , U^) T ^ in R 3 such that ^ . UjQi = and then, for any position (r, a) 
and any N > 1, we would have Ojv(r — v,q)QnU = which contradicts the assumption that the 
solid is not listed in Theorem 4.1. Conversely, if the sequences Q % are R-linearly independent, then 
there exists iVo > 3 such that for all N > Nq, Qm is of rank 3 and the proof is completed. □ 

Proposition 6.4. Assume that the shape Sq of the solid, described by the conformal mapping (3.1), 
is such that e tn ' 2 So = So (the shape of the solid is invariant by rotation of angle ir/2 and center 
0). Then, from the holomorphic potential £ , we can always deduce the values ofr, woe la (the linear 
velocity expressed in a reference fixed frame) and |w| (the absolute value of the rotational velocity). 

Proof. The assumption on the shape Sq means that /(H) = f(£l) where / is the conformal mapping 
defined by: 



with Cfc = e iw ' 2 Ck for any k 
we obtain that: 



f(z) :^c lZ + ^c k z k , (zefl), 
fc<i 

1 and k < — 1. Replacing the function / by / in the computations, 



\i(v) 



E 

fc=i 



Aha 



(r- 



A k e^ k+1 ^/ 2 w Q + B k w a e^ k - 1 ^' 2 + noC k e lk ^ 2 



(6.7) 



and the coefficients X k (v), defined equivalently by (6.2) and by (6.7), must be equal for all r, v £ C, 
a G R/27T and all wq G C and u G R. In particular, for r = v = and a = 0, we obtain 
that A n {\ - e l( * n+1 ^/ 2 )w - S n (l - e^-^^wo = and C„(l - e inn / 2 )uj = for all w G C, 
u € R and n > 1. We deduce that (_4„ 7^ 0) «*■ (n = -1 [4]), (B n ^ 0) <&■ (n = 1 [4]) and 
(C n ^ 0)^(7i = [4]). 

Let us consider the problem of detection now we have some extra information about A n , B n 
and C n . Equation (6.6) with n — 2, gives \\{v — r) + A 2 (^) = for all v G C. Two cases have to 
be considered: 

• Either Ai = 0, which means also that \i(y) = for all v € C. It entails, according to 
equation (6.6) with n = 1, that Ki^o = 0. But B\ — \ci\ 2 ^ and hence w — 0. We 
get u y^ otherwise we would have ^ = 0. Now let m be the smallest index such that 
C m 7^ 0. We know that such an index exists, otherwise it would mean that C k = for 
all k > 1 and hence that any rotational motion of the solid generates a complex potential 
equal to 0. This is impossible for a solid not listed in Theorem 4.1. By induction on n 
with equation (6.6), we must have now X n (v) = for all n < m. We use then equation 
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(6.6) with n = m to obtain iu>C m — e~ lma \ m {v) and next equation (6.2) with n = in + 1 
to get \ m +\{v) = me lma (r — i/)iuiC m . Combining these two identities, we eventually obtain 
A m +i(i / ) = m\ m (v)(r — v), hence we deduce the value of r. 



• Or \ 1 ^ 0. In this case, we deduce easily the value of r from the relation \i(v — r) + A 2 (^) = 0. 

Then, equation (6.6) with n = 1 gives us the value of woe la . There exists at least one m > such 
that C m 7^ 0. Once again, equation (6.6) with n = m allows us to deduce the value of \ui\. D 

If we have more information on <So, we can go further in our study: 

Proposition 6.5. In the preceding proposition, assume furthermore that one of the following as- 
sumption is satisfied: 

1. wq 7^ and there exists an integer m such that either A m and A m +i or B m and yB m +4 are 
both different from 0; 

2. uj ^ and there exists an integer m such that C m and C m+ 4 are both different from 0; 

Then, the solid is detectable. 

Proof. We know that we can compute r and WQe' ,a . Assume for instance that the second assumption 
holds (it would be the same reasoning with the first one). In this case, with equation (6.6) specifying 
n = rn and n — m + 4, we can compute uie" na and cje^ m+4 ' Q . We next deduce e l4a and hence the 
orientation of the solid (because of the symmetry property, e lAa suffices to provide the orientation). 
Since m = [4] (because C m ^ 0) we next deduce the values of e lma and uj. We know that there 
exists at least one index p such that A p ^ 0. We use equation (6.6) with n = p to deduce the value 
of A p Woe' lpa or equivalently A p Woe~ wa . Since necessarily p = — 1 [4] then — p = 1 [4] and we can 
deduce the values of both to and e lAa with woe la . □ 

To conclude this section, we give examples of solids without any symmetry, which are detectable 
(one is pictured in Figure 7): 



6.3 Example of detection 

We consider a shape So described by a complex sequence c such that c\ ^ 0, c_4 =^ and c_7 ^ 0, 
all the other coefficients Ck being null. Direct computations lead to: 



k 


1 


2 


3 


4 


5 


6 


7 


8 


A k 











4 
CfC_ 4 








c[c-7 





B k 


Icil 2 














|ci| 2 C 4 C_ 4 








c k 








C^C_4C_7 





2| Cl |^ 4 c_ 4 








(2| Cl | 2 +3|c_ 4 | 2 )c{c_ 7 



Substiting ing these values into System (6.6), we obtain that: 

B lWo = e~ ia Ai, 

= e- l2Q [A 1 ( J /-r) + A 2 H], 



C%iui = e 



-A4W0 = e 



-?3q 



,fc=l v 

Effcli)'"-'-'""'^"" 



(6.8a) 
(6.8b) 

(6.8c) 
(6.8d) 
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• If uj 7^ and w ^ 0: From equation (6.8a) we deduce that Ai 7^ and from (6.8b) we deduce 
the value of r. Equation (6.8c) allows us to determine |w| and 3a up to n. Combining next 
equation (6.8a) and (6.8d), we get 5a. Using Bezout's identity: a = u3a + w5a with u = 2 
and v = — 1, we get a. We determine u>o with equation (6.8a) and ui with equation (6.8c). 

• If u> = 0, wo 7^ 0: We compute r as in the preceding case. Then, we need further calculations: 



B 6 wo = e 



-A'jWq = e 



-ilea 



,fc=l 

7 



fc-1 



E( n , !)<"-' ■.'" '■'Vih 



E(" !)(---r-' £ A feM 



,fc=i 



fc-i 



(6.9a) 
(6.9b) 



With (6.9a) and (6.8a) we get 5a and with (6.9b) and (6.8a) we get 8a. Since 5 and 8 are 
coprime numbers, we deduce the value of a. We conclude as in the preceding case. 

• If uj 7^ and w — 0: With (6.8a) and (6.8b) we deduce that Ai = A 2 (f ) = for all i/eC. 
We rewrite (6.8c) and (6.8d) as 

C 3 iuj = e~"*A 3 , 

= X A {v) + 3(i/ - r)A 3] 

and we deduce first that A3 7^ and then the value of r. We next add the equations: 

5 



Csiui = e 



-i5a 



C s iuj = e- lSa 



£(i: 1 >- r) "~' A > M 

Ufe=l v 7 



(6.10a) 
(6.10b) 



Equations (6.8c) and (6.10a) give us 2a and (6.10a) and (6.10b) give us 3a. Since 2 and 5 
are coprime numbers we get a and then u> with (6.8c). 



7 Tracking 

In this Section, we perform the proof of Theorem 1.4. So we assume that we know the complex 
potential for all t in a time interval [0, T] (T > 0). We have the expression 



£{t,z) 



E^4> i'- -"i "'/■")■ 



i> 



!(*-")'" 



(7.1) 



where Aj(t, 1/) are complex numbers and R{t,v) := limsup .,-_> +00 |Aj(t, i/)! 1 ". At any time i, the 
series is uniformly convergent on {z £ C : |z — v\ > R(t, u)}. For all N > 1, we denote Ajv(t, ^) := 
(Ai(i, z/), . . . , Ajv(£, i/)) T and we have, according to the results of the preceding section: 

G N U(t) = Q N (r(t) - u, («(())-%((, v), 
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Figure 7: Example of a detectable solid as described in Subsection 6.3. 

where we recall that U(t) :— (di(wo(t)),^s(wo(t)),uj(t)) T . In the proof of Proposition 6.3 we have 
shown that if the solid is not one of those described in Theorem 4.1, then there exists N > 1 such 
that Gn has rank 3. It means that there exists (at least) one inverse Q^ allowing one to express 
the velocity as: 

U(t) := g N 1 e N (r(t) - u, ait^Aifit, v). 

We next get: 

w ) = ( o 1 ) Qn @ N(r(t)-v,a(t)) K N {t,v). 

This equation can be rewritten as: 

dt{a(t)) = { iJ^ e ^(K*) -",«(*)) A w (t >V ), 

to which we can apply the Cauchy-Lipschitz Theorem. The proof is then completed. 

8 Conclusion 

In this article, we have proved that not all solids moving in a perfect fluid can be detected by 
measuring the potential of the fluid. This observation has led us to define the notion of detectable 
solids, which is a purely geometric property. When the geometry is described by means of a 
conformal mapping, we were able to exhibit examples of detectable (or partially detectable) solids. 
However, the complete characterization of such solids in terms of the complex sequence (ck)kez 
remains to be done. 
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